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ABSTRACT 

In this work we have investigated the relevant trend of the bolometric correction (BC) at 
the cool-temperature regime of red giant stars and its possible dependence on stellar metallic- 
ity. Our analysis relies on a wide sample of optical-infrared spectroscopic observations, along 
the 3500 A => 2.5/im wavelength range, for a grid of 92 red giant stars in five (3 globular 
+ 2 open) Galactic clusters, along the full metallicity range covered by the bulk of the stars, 
-2.2 [Fe/H] < +0.4. 

Synthetic BVR C I C JHK photometry from the derived spectral energy distributions al- 
lowed us to obtain robust temperature (T c ff) estimates for each star, within ±100 K or less. 
According to the appropriate temperature estimate, black-body extrapolation of the observed 
spectral energy distribution (SED) allowed us to assess the unsampled flux beyond the wave- 
length limits of our survey. For the bulk of our red giants, this fraction amounted to 15% 
of the total bolometric luminosity, a figure that raises up to 30% for the coolest targets 
(Teff 3500 K). Allover, we trust to infer stellar Mboi values with an internal accuracy 
of a few percent. Even neglecting any correction for lost luminosity etc. we would be overes- 
timating Mboi by < 0.3 mag, in the worst cases. Making use of our new database, we provide 
a set of fitting functions for the V and K BC vs. T c ff and vs. (B — V) and (V — K) broad-band 
colors, valid over the interval 3300 K ^ T e ff ^ 5000 K, especially suited for Red Giants. 

The analysis of the BCy and BC^ estimates along the wide range of metallicity spanned 
by our stellar sample show no evident drift with [Fe/H]. Things may be different for the B- 
band correction, where the blanketing effects are more and more severe. A drift of A(B — V) 
vs. [Fe/H] is in fact clearly evident from our data, with metal-poor stars displaying a "bluer" 
(B — V) with respect to the metal -rich sample, for fixed T c g . 

Our empirical bolometric corrections are in good overall agreement with most of the 
existing theoretical and observational determinations, supporting the conclusion that (a) BCk 
from the most recent studies are reliable within < ±0.1 over the whole color/temperature 
range considered in this paper, and (b) the same conclusion apply to BCv only for stars 
warmer than ~ 3800 K. At cooler temperatures the agreement is less general, and MARCS 
models are the only ones providing a satisfactory match to observations, in particular in the 
BC V vs. (B - V) plane. 

Key words: Stars: late-type - Stars: atmospheres - Galaxy: globular clusters: general - 
Galaxy: stellar content - infrared: stars 



* Based on observations made at La Palma, at the Spanish Observatorio del 
Roque de los Muchachos of the IAC, with the Italian Telescopio Nazionale 
Galileo (TNG) operated by the Fundacion Galileo Galilei of INAF. 



1 INTRODUCTION 

A physical assessment of the bolometric emission of stars is a 
mandatory step for any attempt to self-consistently link observa- 
tions and theoretical predictions of stellar evolution. The impor- 
tance of this comparison actually reverberates into a wide range of 
primary astrophysical questions, ranging from the validation of the 
reference input physics for nuclear reactions in the stellar interiors 



2 A. Buzzoni et al. 



to the study of integrated spectrophotometric properties of distant 
galaxies, through stellar population synthesis models. 

By definition, the effective temperature (Toff) and physical 
size (R) of a star provide the natural constraint to its emerging 
flux, as L oc R 2 T* s . If L is a known property for a star, then we 
could physically "rescale" the spectral energy distribution (SED), 
and infer, from the observed flux, the distance of the body, d, or its 
absolute size (R), through a meas ure of the apparent angular exten- 



sion, 6 = (R/df oc LTrJ d~ 2 dRidgwav et alJ|l98d ; lDvck et al.1 
ll996j ; |Perrin et alll998l ; lRichichi et alj|l998l) . 

As well known, however, L cannot, in principle, be directly 
measured, requiring for this task an ideal detector equally sensitive 
to the whole spectral range. The lack of this crucial piece of infor- 
mation is often palliated by indirect observing methods, trying to 
pick up the bulk of stellar emission through broad-band photometry 
within the appropriate spectral r ange according to target tempera- 
ture^ Relying on this approach, Johnson (1966) deriv ed the bolo- 
metric vs. temperature scale for red giant stars, while ICode et all 
dl976h explored the same relation for hot early-type stars, through 
satellite-borne UV observations. As an alternative way, many au- 
thors tried a fully theoretical assessment of the problem, by study- 
ing the /boi vs. fx relationship on the basis of model grids of stellar 
atmospheres and replacing observations with synthetic photome- 
try directly computed on the t heoretical SED dBertone et alj|2004l ; 
IBessell, Castelli. & Plezll 19981) . 

Rather than focussing on luminosity, Wesselink ( 1969) 
originally proposed a further application of this method, 
just looking at the bolometric surface brightness, namely 
/i = fboi/8 2 , to lead to a refined temperature scale of stars 
in force of the fundamental relationship /i = a (a 
being the Stefan-Boltzmann constant). The so-called surface- 
brightness technique, then better recognized as the IR-flux 
method (IRFM), has been ex tensively applied to the study of 



red giant and supergian t stars (fBlackwell, Shallis, & Selbv 
di Be nedetto & Rabbial fl987l; iBlackwell & Lvnas-Gravl 



1979 



1994; 



Alonso, Arri bas. & Marti'nez-Roger 1999; | Rarmrez & Melendezl 
20051 ; iGonzalez Hernandez & Bonifacio 2009) taking advantage 



of its distance-independent results, providing to match the angular 
measure of stellar radii with the estimate of the bolometric flux 
from infrared observations, i.e. n — (/boi//m)/iR- 

Although in different forms, all the previous methods used 
theoretical models of stellar atmospheres to derive the appropriate 
"correcting factor" TZ = /boi//(A) and convert observed or syn- 
thetic monochromatic magnitudes, m(A) to the bolometric scale0 
Taking the Sun as a reference source for our calibration, we could 
write more explicitely: 



[rn hol - m(A)] - [m bol - m(A)] = -2.5 log(7£/7£„) 



(1) 



Equation Q} actually leads to the straight definition of bolo- 
metric correction, BC(X), namely 



BC{\) = (m bo i - m(A)) = -2.5 \oglZ + BC(X)@ 



(2) 



Aside from the historical definition, that originally con sidered 
BC only to photographic (m pg ) or visual (my) magnitudes ( Kuiper 



1 Recalling that emission peak roughly obeys the Wien law, i.e. T\ pc ^k — 
const. 

2 To a more detailed analysis, note that the ratio TZ dimensionally matches 
the definition of "equivalent width", and it gives a measure of how "broad" 
is the whole SED compared to the monochromatic emission density at the 
reference A. 



Il938h . one can nowadays easily extend the definition to any wave- 
band. A careful analysis of eq. ([2} makes clear some important 
properties of BC: i) the value of TZ is a composite function of stel- 
lar fundamental parameters, namely TZ = TZ(T c a, log g, [X/H]) 
so that, for fixed effective temperature, BC may display some de- 
pendence on stellar gravity (g) and chemical composition ([X/H]); 
ii) the value of TZ (and, accordingly, of BC) is minimum when our 
observations catch the bulk of stellar luminosity. For this reason, 
high values of BCv must be expected when observing for instance 
cool giant stars in the V band, or hot O-B stars in the infrared K 
band. Hi) The definition of the BC scale strictly depends on the as- 
sumed reference value for the S un, that therefore must univoc ally 
fix the "zero point" of the scale ( iBessell. Castelli, & Pledl998l) . 

In this framework, we want to tackle here the central ques- 
tion of the possible BC dependence on stellar metallicity. This ef- 
fect could be of special importance, in fact, in order to more confi- 
dently set the bolometric vs. temperature scale for cool red giants, 
where the intervening absorption of diatomic (TiO in primis) and 
triatomic (H2O) molecules heavily modulate the stell ar SED with 
sizeable effects on optical and NIR magnitudes (e.g. lGratton et al.l 
Il982l ; iBertone etid] I2OO8). As a matter of fact, still nowadays 
the many efforts devoted to the definition of the BC vs. log T c n 
relationship led to non univocal conclusions, with large discrep- 
ancies among the different sources in the litera ture as far as 
spectral type or later are concerned ^Flowed 
ll & Wood 1 1984 iHoudashelt. Bell, & Sweigartl 



stars of K 



1975 



2000 



1977; Bessell 

Bertone et al.ll2004l ; IWorthev & Leel2006l) . 

This issue has actually an even more important impact on the 
study of the integrated spectrophotometric properties of resolved 
and unresolved stellar systems, as red giants and other Post-main 
sequence (P MS) stars provide a prevailing fraction (2/3 or more, 
iBuzzonil 19891) of the total luminosity of the population. A fair defi- 
nition of the BC scale becomes therefore of paramount importance 
to self-consistently convert the oretical H-R and observed c-m dia- 
gram s of a stellar population dFloweiiri996l ; IVan denBerg & Cleml 
2003) and more confidently assess the physical contribution of the 
different stellar classes. 

A study of the BC dependence on metal abundance has been 
previously attempted by many authors mainly relying on a fully 
theoretical point of view to exploit the obvious advantage of stel- 
lar models to account in a controlled way for a global or selective 
chang e of metal abunda nce. In this regard. jTripicco & Belli dl995l) 
and lCassisi et al.l d2004l) . among others, tried to explore the effect 
of a eleme nts enhancement (na mely O, Mg, Ca, Ti etc.) in stellar 
SED, while lGirardi et al.l {2007) focussed on the possible impact of 
Helium abundance on BC. As a major drawback of these efforts, 
however, one has to report the admitted limit of model atmospheres 
in accurately describe the spectrophotometric properties of K- and 
M-type stars, that are cooler than 4000 K (see lBertone et alJ[20o3 . 
on this important point). 

On the other hand, a fully empi rical approach 

has been devised by [M ontegriffo et alj (l998) and 
lAlonso. Arribas. & Martmez-RogeJ dl999l) , among others, 
trying to reconstruct stellar SED, and therefrom infer the bolo- 
metric flux, /boi, through optical broad-band photometry of stars 
in the Galactic field or in globular clusters. A recognized limit of 
these studies is, however, that they may suffer from the lack of 
coverage of the stellar parameter space offered by the observations. 
Moreover, as far as the cool-star sequence is concerned, optical 
multicolor photometry, alone, partially misses the bulk of stellar 
emission (more centered toward the NIR spectral window); in ad- 
dition, by converting broad-band magnitudes into monochromatic 



flux densities, the stellar SED is reconstructed at very poor spectral 
resolution, thus possibly loosing important features that may bias 
the inferred bolometric energy budget. 

On this line, however, we want to further improve the analy- 
sis proposing here more complete spectroscopic observations for 
a large grid of red-giant stars in several Galactic clusters along 
the entire metallicity scale from very metal-poor (i.e. [Fe/H] ~ 
—2.2 dex) to super-solar ([Fe/H] ~ +0.4 dex) stellar popula- 
tions. Our observations span the whole optical and NIR wavelength 
range, thus allowing a quite accurate shaping of stellar SED. As we 
will demonstrate in the following of our discussion, our procedure 
allowed us to sample about 70-90 % of the total emission of our 
sample stars, thus leading to a virtually direct measure of /bob even 
for M-type stars as cool as 3500 K. 

We will arrange our discussion by presenting, in Sec. 2, our 
stellar database together with further available information in the 
literature. The analysis of the observing material will be assessed 
in more detail in Sec. 3, while in Sec. 4 we will derive the SED for 
the whole sample leading to an estimate of the effective tempera- 
ture and bolometric correction for each star. The discussion of the 
inferred BC-color-temperature scale will be the focus of Sec. 5, es- 
pecially addressing the possible dependence of BC on stellar metal- 
licity. The comparison of our results with other relevant BC calibra- 
tion in the literature will also be carried out in this section, while in 
Sec. 6 we will summarize the main conclusions of our work. 
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2 CLUSTER DATABASE SELECTION 

As we mainly aim at probing the impact of metallicity on the BC 
of stars at the low-temperature regime, a demanding constraint to 
set up our target sample was to explore a range as wide as possible 
in [Fe/H], and pick up red giant stars with accurate measurements 
of their metallicity. The cluster population in the Galaxy naturally 
provided the ideal environment for our task. By combining globular 
and open clusters one can easily span the whole metallicity range 
pertinent to Pop I and II stars in our and in external galaxies. We 
therefore selected five template systems, namely the three metal- 
poor globular clusters M15, M2 and M71, and two metal-rich open 
clusters NGC 188 and NGC 6791 such as to let metallicity span 
almost three orders of magnitude, from [Fe/H] = —2.3 up to +0.4. 

For each cluster, a subset of ~ 20 suitable targets have then 
been identified among the bright est and coolest red g iants from 
the 2MASS infrared c-m diagram dSkrutskie et alj2006h . In assem- 
bling the dataset we also took care of picking up those objects out 
of more severely crowded regions of the clusters, and clearly recog- 
nizable in bright asterisms such as to reduce the chance of misiden- 
tification at the telescope. 

The final set of target stars is summarized, for each cluster, in 
the five panels of Fig.[T]and in the series of TablesQ]to[5] We eventu- 
ally considered 92 stars in total, of which 21 are in M15, 18 in M2, 
17 in M71, 16 in NGC 188, and 20 in NGC 6791, respectively. For 
each star, the tables always report the 2MASS id number (col. 1) 
and the alternative cross-identification, according to other reference 
photometric catalogs, when available. The 2MASS J2000 coordi- 
nates on the sky and the corresponding J, H, K magnitudes are 
also always reported, together with a compilation of B,V, R c , I c 
observed magnitudes according to the best reference catalogs for 
each cluster, as reported in the literature. When required, dered- 
dened apparent magnitudes have been computed according to the 
color excess E( B — V) as labelled in the header of each table. 



Figure 1. Apparent c-m diagram of the five clusters included in our analysis, 
according to 2MASS J and K photometry. Big squares along the red giant 
branch mark the selected targets in our sample. 

3 OBSERVATIONS AND DATA REDUCTION 

Spectroscopic observations of our stellar sample have been col- 
lected during several runs between June and October 2003 at the 
3.5m Telecopio Nazionale Galileo (TNG) of the Roque de los 
Muchachos Observatory, at La Palma (Canary Islands, Spain). A 
summary of the logbook can be found in Table|6] 

Optical spectroscopy was carried out with the LRS FOSC 
camera; a composite spectrum was collected for each target by 
matching a blue (grism LRB along the AA3500 — 8800 A wave- 
length rangefj] and a red setup (grism LRR, between AA4500 — 
10300 A). In both cases the grisms provided a dispersion of 
2.8 A px _1 on a 2048 x 2048 thinned and back-illuminated Loral 
CCD, with a 13.5/im pixel size. In order to collect the entire flux 
from target stars, we observed through a 5"wide slit; this condition 
actually made spectral resolution to be eventually constrained by 
the seeing figure (typically about 1-1. 5" along the different nights), 
thus ranging between 10 and 15 A (FWHM). This is equivalent 
to a value of R = A/ A A between 600-1000. Whenever possible, 
and avoiding severe crowding conditions of the target fields, the 
longslit was located at the parallactic angle. Wavelength calibration 
and data reductions were performed following standard procedures. 

The optical spectra have then been accompanied by the corre- 
sponding observations taken at infrared wavelength with the NICS 
camera at the Nasmyth focus of the TNG. The camera was coupled 
with a Rockwell 1024 x 1024 Hawaii- 1 HgCdTe detector. We took 
advantage of NICS unique design using the Amici grism coupled 

3 Although nominally extended to 9500 A, LRB spectra result severely af- 
fected by second-order spectral emission in their red tail. For this reason, 
during data reduction, spectra have been clipped retaining only the wave- 
length region blueward of 8800 A. 
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Table 1. Cluster properties and stellar database for cluster M 15 



M 15: 



E(B-V) = 0.10 



[Fe/H] = -2.26 



ID 








a 


S 


B 




V 




Ic 


J 


H 


K 


(a) 


to 


(c) 




(J2000.0) 


to 


(6) 


to 


to 


to 


(a) 


(a) 


(a) 


21300002+1209182 


165 


71 


21 


30:00.02 


12:09:18.24 


15.334 


14.395 


14.3460 


13.330 


13.2709 


12.479 


11.926 


11.824 


21295705+1208531 


959 


6 


21 


29:57.06 


12:08:53.11 


14.549 


13.426 


13.4946 


12.144 


12.2129 


11.282 


10.691 


10.573 


21295532+1210327 


337 


60 


21 


29:55.33 


12:10:32.80 


15.229 


14.313 


14.3694 


13.165 


13.2132 


12.452 


11.899 


11.786 


z i _5uuuyu+ j zuo j 1 1 








30:00.91 


12:08:57.13 


















21295473+1208592 


330 


25 


21 


29:54.73 


12:08:59.24 


14.821 


13.691 


13.7581 


12.444 


12.5006 


11.591 


11.065 


10.906 


21300461+1210327 


369 




21 


30:04.62 


12:10:32.73 


14.851 


13.836 








11.858 


11.272 


11.165 


21295560+1212422 


533 


665 


21 


29:55.61 


12:12:42.29 


14.562 


13.459 


13.5218 




12.2237 


11.336 


10.723 


10.609 


21300514+1210041 


372 




21 


30:05.15 


12:10:04.18 


15.186 


14.288 








12.430 


11.929 


11.776 


21295836+1209020 




166 


21 


29:58.37 


12:09:02.01 






13.8205 




12.5987 


11.700 


11.112 


11.042 


21295618+1210179 




631 


21 


29:56.18 


12:10:17.93 






12.7694 




11.3768 


10.414 


9.781 


9.649 


21295739+1209056 




7 


21 


29:57.39 


12:09:05.69 






13.7397 




12.5054 


11.632 


11.070 


10.948 


21300097+1210375 




65 


21 


30:00.98 


12:10:37.60 






13.8739 




12.6289 


11.726 


11.171 


11.017 


21300431+1210561 


368 




21 


30:04.32 


12:10:56.16 


14.649 


13.559 








11.459 


10.893 


10.757 


21301049+1210061 


621 




21 


30:10.49 


12:10:06.18 


14.563 


13.406 








11.151 


10.562 


10.438 


21300739+1210330 


604 




21 


30:07.40 


12:10:33.06 


14.961 


13.986 








11.964 


11.399 


11.264 


21300569+1210156 






21 


30:05.70 


12:10:15.68 












12.156 


11.596 


11.480 


21300553+1208553 






21 


30:05.54 


12:08:55.35 












12.357 


11.835 


11.719 


21295756+1209438 






21 


29:57.57 


12:09:43.85 












10.096 


9.429 


9.269 


21295082+1211301 






21 


29:50.83 


12:11:30.18 












11.326 


10.725 


10.612 


21295881+1209285 




59 


21 


29:58.82 


12:09:28.59 






14.5465 




13.5061 


11.088 


10.568 


10.353 


21295716+1209175 




273 


21 


29:57.17 


12:09:17.52 






13.1662 




11.7880 


10.867 


10.220 


10.112 



( a > fro m 2MASS; 

< b > from Cohen, Briley, & St etsot] feOOSh : 
( c ' from 



Rosen berg etalJfeOOOh 



Table 2. Cluster properties and stellar database for cluster M 2( a ' 



M2: 

ID 



E(B-V) = 0.06 



[Fe/H] = -1.62 



(J2000.0) 



21333827-0054569 


21 


33:38.28 


-00:54:56.92 


10 


542 


9 


827 


9.672 


21333095-0052154 


21 


33:30.96 


-00:52:15.47 


1 1 


568 


10 


952 


10.814 


21332468-0044252 


21 


33:24.69 


-00:44:25.21 


12 


549 


12 


006 


11.886 


21331771-0047273 


21 


33:17.71 


-00:47:27.31 


10 


(.65 


9 


961 


9.821 


21331723-0048171 


21 


33:17.24 


-00:48:17.10 


1 1 


1 12 


10 


429 


10.301 


21331790-0048198 


21 


33:17.91 


-00:48:19.82 


1 1 


746 


1 1 


103 


11.017 


21331854-0051563 


21 


33:18.55 


-00:51:56.33 


1 1 


779 


1 1 


137 


11.019 


21331948-0051034 


21 


33:19.49 


-00:51:03.42 


1 1 


963 


1 1 


299 


11.214 


21331923-0049058 


21 


33:19.23 


-00:49:05.84 


12 


2S0 


1 1 


695 


11.579 


21332588-0046004 


21 


33:25.89 


-00:46:00.44 


12 


313 


1 1 


756 


11.600 


21333668-0051058 


21 


33:36.68 


-00:51:05.89 


10 


730 


10 


026 


9.880 


21333520-0046089 


21 


33:35.21 


-00:46:08.91 


10 


993 


10 


324 


10.174 


21333488-0047572 


21 


33:34.88 


-00:47:57.25 


1 1 


265 


10 


589 


10.455 


21333593-0049224 


21 


33:35.94 


-00:49:22.44 


1 1 


420 


10 


750 


10.650 


21333432-0051285 


21 


33:34.33 


-00:51:28.50 


1 1 


490 


10 


828 


10.722 


21332531-0052511 


21 


33:25.32 


-00:52:51.17 


1 1 


938 


1 1 


300 


11.203 


21333109-0054522 


21 


33:31.09 


-00:54:52.28 


12 


086 


1 1 


526 


11.376 


21333507-0051097 


21 


33:35.07 


-00:51:09.72 


12 


609 


12 


056 


11.962 



(°) all the data are from 2MASS; 



with two slits 0.5 "and 5 "wide, the latter being used for a complete 
flux sampling of the target stars. The spectra cover the entire wave- 
length range from 8000 A to 2.5 fim at a resolving power (for a 
0.5"slit) which varies between R — 80 and 140 along the spec- 
trum. In acquiring spectra, background subtraction and flat-fielding 
correction were eased by a standard dithering procedure on tar- 
get images, while the wavelength calibration directly derived from 
the standard reference table providing the dispersion relation of the 
system. The MIDAS ESO package, and specifically its LONGSLIT 
routine set, has been used for the whole reduction procedure, both 
for optical and infrared spectra. 



3.1 Flux calibration 

Given the nature of our investigation, special care has been de- 
voted to suitably fluxing both optical and infrared spectra. This 
has been carried out by repeated observations, both with LRS and 
NICS, of a grid of sp ectr ophotometric stand ard stars from the list of 
iMassev et all < ll988» and lkunt et all dl998h . as reported in Table[6] 
Note, however, that the lack of an appropriate SED calibration of 
standard stars along the entire wavelength range of our observations 
required a two-step procedure, relying on the direct observation 
of Veg a as a primary calibrator, according to Tokunaga & Vaccal 
d2005» results Given the outstanding luminosity of this star we had 
to observe through a 10 mag neutral filter to avoid CCD saturation, 
and create a secondary calibrator (namely HD192281) observed 
both with and without the neutral density filter. 

Concerning the applied correction for atmosphere absorption, 
we had to manage two delicate problems. From one hand, in fact, 
the intervening action of Sahara dust (the so-called "calima effect") 
may abruptly increase the atmosphere opacity at optical wave- 
length. This is a recurrent feature for summer nights at La Palma, 
and it can severily affect the observing output, especially when 
dealing with absolute flux calibration. A careful check with re- 
peated observations of the same standard stars along each night 
allowed us to assess the presence of dust in the air. This confirmed, 
for instance, that along our observing runs, the night of Aug 07, 
2003 displayed an outstanding (i.e. a factor of four higher than the 
average) dust extinction. 

On the other hand, atmosphere water vapour can also play 
a role by affecting in unpredictable ways the infrared obser- 
vations. Telluric H2O bands about 1.10, 1.38 and 1.88 ^tm 
dFuensalida & Alonsol ll998; Manduca & Belllll979l) . just restrain- 
ing to the Amici wavelength range, may in fact strongly contami- 
nate the intrinsic H2O absor ption bands of stella r SED, especially 
for stars cooler than 3500 K ( iBertone et aT1 l2008). This effect may 
act on short timescales along the night, so that it cannot be recon- 
ducted to an average nightly extinction curve, as for optical obser- 
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Table 3. Cluster properties and stellar database for cluster M 7 1 



M 71: E(B-V) = 0.25 [Fe/H] = -0.73 



ID 






a 


5 


B 




V 




1 


H 


K 


(a) 


(M 


(c) 


(J2000.0) 


00 


(6) 


O) 


(c) 


(a) 


(a) 


(a) 


19535325+1846471 


2672 


256 


19:53:53.25 


18:46:47.13 


13.905 


12.314 


12.2085 


10.3988 


9.090 


8.197 


8.040 


19534750+1846169 


2222 


540 


19:53:47.51 


18:46:16.99 


14.431 


13.137 


13.0010 


11.5156 


10.452 


9.698 


9.588 


19535150+1848059 


2541 


892 


19:53:51.50 


18:48:05.91 


14.079 


12.436 


12.3250 


10.4275 


9.079 


8.207 


7.968 


19535064+1849075 


2461 


331 


19:53:50.64 


18:49:07.52 


14.466 


13.064 


12.9955 


11.4204 


10.215 


9.446 


9.271 


19534575+1847547 


2079 


648 


19:53:45.76 


18:47:54.80 


14.247 


12.606 


12.4924 


10.5109 


9.094 


8.203 


7.974 


19534827+1848021 


2281 


309 


19:53:48.27 


18:48:02.17 


14.078 


12.492 


12.3636 


10.5500 


9.177 


8.270 


8.094 


19534656+1847441 


2145 


46 


19:53:46.57 


18:47:44.19 


14.838 


13.623 


13.5524 


12.2176 


11.228 


10.569 


10.435 


19535369+1846039 


271 1 


172 


19:53:53.70 


18:46:03.98 


15.527 


14.578 


14.4974 


13.3402 


12.500 


11.998 


11.896 


19534905+1846003 


2337 


303 


19:53:49.05 


18:46:00.34 


14.601 


13.410 


13.3436 


11.9991 


10.950 


10.276 


10.186 


19534916+1846512 


2347 


6 


19:53:49.16 


18:46:51.22 


14.997 


13.709 


13.6219 


12.2031 


11.151 


10.434 


10.301 


19534178+1848384 


1772 




19:53:41.79 


18:48:38.46 


15.877 


14.694 






12.183 


11.521 


11.402 


19535676+1845399 


2921 




19:53:56.77 


18:45:39.95 


15.747 


14.605 






12.197 


11.529 


11.455 


19533962+1848569 


161 1 




19:53:39.62 


18:48:56.99 


15.695 


14.627 






12.494 


11.974 


11.888 


19533864+1847554 


1543 




19:53:38.64 


18:47:55.45 


15.475 


14.222 






11.751 


11.151 


11.037 


19534615+1847261 




580 


19:53:46.15 


18:47:26.11 






13.1140 


11.5109 


10.336 


9.543 


9.395 


19535610+1847167 


2885 


1066 


19:53:56.10 


18:47:16.76 


13.577 


11.905 


12.4009 


9.2167 


7.943 


7.078 


6.681 


19534941+1844269 


2365 




19:53:49.41 


18:44:26.98 


13.863 


12.107 






8.058 


7.105 


6.863 



W fro m2MASS; 

( 6 > from Geffert & Maintzl feOOOl) ; 
( c ' from Rosenber g et a l. 1 2000) 



Table 4. Cluster properties and stellar database for NGC 188 



NGC 188: E(B-V) = 0.082 [Fe/H] = -0.02 

ID a S B V R c I c J H K 

(a) (b) (c) (J2000.0) W < c > < b > < c > < c > ( c > W (") M 



00445253+8514055 


466S 


N18 


3-1-69 


00:44:52.54 


85 


14:05.54 


13.613 


13.579 


12.319 


12.357 


11.598 


11.087 


10.098 


9.461 


9.304 


00475922+8511322 


5887 


N18 


3-II-181 


00:47:59.23 


85 


11:32.28 


13.587 


13.428 


12.135 


12.197 


11.429 


10.894 


9.891 


9.203 


9.100 


00465966+8513157 


5085 


N18 


8-1-105 


00:46:59.66 


85 


13:15.71 


13.603 


13.538 


12.362 


12.422 


11.732 


11.269 


10.349 


9.789 


9.639 


00453697+8515084 


5927 


N18 


3-1-57 


00:45:36.97 


85 


15:08.43 


14.799 


14.760 


13.658 


13.706 


13.039 


12.571 


11.709 


11.149 


11.024 


00442946+8515093 


4636 


N18 


3-1-59 


00:44:29.46 


85 


15:09.39 


14.986 


14.950 


14.005 


14.046 


13.385 


12.962 


12.202 


11.653 


11.520 


00473222+8511024 


5133 


N18 


3-II-187 


00:47:32.22 


85 


11:02.45 


15.171 


15.132 


14.077 


14.140 


13.490 




12.234 


11.700 


11.567 


00554526+8512209 


6175 






00:55:45.27 


85 


12:20.92 


12.224 




10.834 








8.441 


7.631 


7.520 


00463920+8523336 


4X43 






00:46:39.21 


85 


23:33.67 


12.890 




11.569 








9.292 


8.597 


8.441 


00472975+8524140 


4829 






00:47:29.76 


85 


24:14.09 


13.965 




12.781 








10.783 


10.210 


10.114 


00441241+8509312 


4756 






00:44:12.42 


85 


09:31.23 


12.933 




11.404 








8.580 


7.892 


7.652 


00432696+8509175 


4408 






00:43:26.96 


85 


09:17.58 


14.242 




13.199 








11.293 


10.706 


10.591 


00471847+8519456 


4909 






00:47:18.48 


85 


19:45.65 


14.255 




13.010 








10.908 


10.289 


10.187 


00461981+8520086 


4524 






00:46:19.81 


85 


20:08.61 


13.663 




12.468 








10.385 


9.816 


9.674 


00463004+8511518 


5894 






00:46:30.05 


85 


11:51.89 


15.142 




14.052 








12.185 


11.695 


11.518 


00490560+8526077 


5835 






00:49:05.60 


85 


26:07.77 


13.921 




12.717 








10.594 


9.956 


9.825 


00420323+8520492 


= SAO 109 




00:42:03.23 


85 


20:49.23 


11.40 


(d) 


9.89 


(d) 






7.064 


6.387 


6.130 



M from 2MASS: 

( & ) froml Platais etai]fe003l) 
( c ) from 



Stetson. McClure. & VandenBerd J2004h ; 



< d ) B and V photometry from SIMBAD 



vations. The H2O contamination in each spectrum was therefore 
corrected by re-scaling the average extinction curve to minimize 
the residual water vapour feature in the stellar spectra. 



Allover, the full calibration procedure led us to consistently 
assemble the LRB-LRR-Amici spectral branches and obtain a nom- 
inal SED of target stars along the 3450-25000 A wavelength range. 
However, just an eye inspection of the full spectra made evident 
in some cases a residual systematic component causing a "glitch" 
at the boundary connection between LRS and NICS observations. 
Clearly, this effect urged us to further refine our analysis taking into 
account the supplementary photometric piece of information, as we 
will discuss in more detail in the next section. 



3.2 Photometry and spectral "fine tuning" 

The relevant database of broad-band photometry available in the 
literature for all stars in our sample can be usefully accounted for 
our analysis as a supplementary tool to tackle the inherent difficulty 
in reproducing the overall shape of stellar SED at the required ac- 
curacy level over the entire range of our observations. 

As summarized in Tables Q] to [5] a wide collection of photo- 
metric catalogs can be considered, providing multicolor photome- 
try along the range spanned by LRS and NICS spectra. Facing the 
observed values, one can similarly derive a corresponding set of 
multicolor synthetic magnitudes relying on the assembled SED of 
each star. Operationally, from our /(A) values we need to numeri- 
cally assess the quantity 
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Table 5. Cluster properties and stellar database for cluster NGC 679 1 



NGC 6791: E(B-V) = 0.117 [Fe/H] = +0.4 



ID 

(a) 


(&) 


(<0 


a 5 
(J2000.0) 


B 

00 


(<=) 


V 

W 


(<0 


lr 

(c) 


J 

(a) 


H 

(a) 


K 

(a) 


19210807+3747494 


6697 


3475 


19:21:08.07 


37:47:49.41 


15.279 


15.275 


13.909 


13.978 


12.668 


11.675 


11.071 


10.919 


19204971+3743426 


10807 


3502 


19:20:49.72 


37:43:42.67 


15.554 


15.563 


13.956 


13.982 


10.990 


9.041 


8.167 


7.815 








19:20:52.60 


37:44:28.18 


















19205580+3742307 


11799 


3574 


19:20:55.81 


37:42:30.75 


16.307 


16.297 


14.934 


14.957 


13.609 


12.622 


11.993 


11.860 


19205671+3743074 


11308 


2478 


19:20:56.72 


37:43:07.46 


16.000 


15.984 


14.633 


14.660 


13.325 


12.351 


11.756 


11.586 


19210112+3742134 


12010 


3407 


19:21:01.12 


37:42:13.45 


15.942 


15.928 


14.433 


14.455 


12.901 


11.821 


11.130 


10.938 


19211606+3746462 


7750 




19:21:16.06 


37:46:46.26 


15.472 




13.871 






8.914 


8.053 


7.714 


19213656+3740376 


12650 




19:21:36.56 


37:40:37.63 


15.727 




14.174 






11.431 


10.635 


10.438 


19210326+3741190 


13637 




19:21:03.27 


37:41:19.04 


15.722 




14.348 






12.120 


11.516 


11.362 


19213635+3739445 


13082 




19:21:36.36 


37:39:44.57 


16.186 




14.825 






12.449 


11.728 


11.608 


19212437+3735402 


15790 




19:21:24.37 


37:35:40.29 


15.837 




14.442 






12.134 


11.546 


11.354 


19212674+3735186 






19:21:26.75 


37:35:18.60 












11.622 


10.925 


10.735 


19211632+3752154 


3254 




19:21:16.32 


37:52:15.46 


15.282 




13.998 






11.776 


11.068 


10.967 


19211176+3752459 


2970 




19:21:11.76 


37:52:46.00 


15.676 




14.336 






12.174 


11.557 


11.420 


19202345+3754578 


1829 




19:20:23.45 


37:54:57.82 


14.592 




12.866 






8.029 


7.133 


6.787 


19205149+3739334 






19:20:51.50 


37:39:33.44 












7.356 


6.516 


6.201 


19203285+3753488 


2394 




19:20:32.85 


37:53:48.87 


15.056 




13.417 






8.463 


7.535 


7.224 


19200641+3744452 


9800 




19:20:06.42 


37:44:45.28 


14.670 




13.307 






10.831 


10.094 


9.943 


19200882+3744317 


10034 




19:20:08.83 


37:44:31.71 


15.353 




13.710 






7.916 


6.989 


6.670 


19203219+3744208 


10223 




19:20:32.20 


37:44:20.81 


16.421 




14.854 






8.176 


7.262 


6.874 



( a > fro m2MASS; 

( 6 ) froml Kaluznv & Rucinskil J 19951): 
( c ' from 



Stetson. Bruntt. & Grundah] 120031) 



Table 6. Logbook of TNG observations along 2003 



Obs. date 
(2003) 


Instrument 


Targets 


Standards ' 


Jul 29 


LRS 


NGC6791 


HD 192281 


Jul 30 


LRS 


NGC6791 


HD192281, SAO48300, WOLF1346 


Jul 31 


LRS 


NGC6791 


HD192281, SAO48300, WOLF1346 


Aug 6 


LRS 


M71 


HD192281, SAO48300, WOLF1346 


Aug 7 


LRS 


MIS 


HD192281, SAO48300, WOLF1346 


Aug 11 


NICS 




HD192281,SAO48300 


Aug 12 


NICS 




Vega 


Aug 18 


NICS 


M7I 


HD192281 


Aug 19 


NICS 




Vega 


Aug 20 


NICS 


M15, M71, 


HD192281, SAO48300, WOLF1346 






NGC188 


Vega 


Aug 21 


LRS 


M2 


HD192281 


Aug 23 


LRS 


M2, M15 


HD 192281 


Aug 26 


LRS 


NGC188 




Aug 27 


LRS 


NGC188 


HD 192281 


Aug 31 


NICS 


M71 


HD192281 


Sep 1 


NICS 


M71.NGC6791 


HD 192281 


Sep 3 


NICS 


M2, M15, NGC188 


SAO48300 


Sep 4 


NICS 


NGC 188 


HD 192281 


Sep 5 


NICS 


NGC188, NGC6791 


HD 192281 


Oct 14 


NICS 


MI5 


HD192281 


Oct 15 


NICS 


M2 


HD192281 



j_HD192281 and WOLF1346 from optical ralihrarin n hv iMassev et al I 
ll98Sl) : S AO4S300 from THK nhntnrnetrir calibration bv lHunt et alJ l 1998): 
Vega from lTokunaga & Va cca 12005) 



tudes are available for most clusters, observed colors have been 
computed for each available V dataset and are displayed with a 
different marker (either dot or square) in the plots. 

Just a glance to Fig. [2] makes evident that systematic offsets 
are present between observed photometry and synthetic magni- 
tudes. This may partly be due to zero-point uncertainty in comput- 
ing eq. 10, as well as to residual systematic drifts inherent to our 
spectral flux calibration. In addition, from the figure one has also 
to report a few outliers in every band, and a notably skewed distri- 
bution of B residuals. To recover for this systematics we devised 
an iterative 3<j clipping procedure on the data of Fig. [2] to reject 
deviant stars and lead synthetic magnitudes to match the standard 
photometric system of the observed catalogs. Our results are dis- 
played in graphical form in the plots of Fig. [3] 

After just a few rejections, our procedure quickly converged 
to mean magnitude offsets ((Obs — Syn), see Table O to correct 
eq. <[3j output. After correction for this systematics, our final syn- 
thetic photometry of cluster stars (not accouting for Galactic red- 
dening) is collected in Tables [8] and [9] According to Table [7] note 
that a a = 0.095 mag in total magnitude residuals evidently im- 
plies an internal accuracy in our spectral flux calibration of target 
stars better than 10%. 



f f(x)s(\y d\ 

being m J syn the synthetic magnitude in the "j-th" photometric 
band, identified by a filter response S(X) J and a calibrating zero- 
point flux fi . For our calculations we relied on the lBuzzonil ( 2005) 
reference data (see Table 1 therein). 

A comparison of our output with the available photometry is 
displayed in Fig. [2] The magnitude difference (in the sense "syn- 
thetic" - "observed"), is plotted in the different panels of the figure 
vs. observed color, according to the different photometric catalogs 
quoted in Tables Q] to [5] As typically two sources for V magni- 



3.3 Stellar outliers 

It could be interesting to analyze in some detail the deviant stars in 
our Am clipping procedure in order to collect further clues about 
their nature. Apart from the obvious impact of photometric errors, 
3<J outliers may in fact more likely be displaying signs of an intrin- 
sic physical variability in their luminosity. 

As summarized in Table 1 101 in total 10 stars have been found 
to significantly (> 3<r) deviate from the literature compilations. A 
careful check of their identifications on the SlMBAD database in- 
dicates that at least 3 of them are known variable (typically semir- 
regulars or irregulars, as expected for their nature of late-type red 
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Figure 2. The B,V, R c , Ic, J, H, K magnitude residuals between syn- 
thetic and observed magnitudes (in the sense "syn" - "obs") for the 94 stars 
in our sample. plotted vs. literature colors, according to the data of Table[T] 
to [5] Synthetic magnitudes derived from the numerical integration of the 
observed SED with the Johnson-Cousins filters. All the available photom- 
etry has been accounted for. Some stars with multiple V datasets appear 
therefore twice in the plots and are singled out by dot and square markers, 
respectively. 
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Figure 3. The histogram of magnitude residuals for the data of Fig.|2] after 
correction for the systematic offsets, according to Table [7] A total of 492 
measures have been accounted for, as labelled in the global histogram of 
the bottom panel, including multiple photometry sources in the literature 
from Tables [T] to \5\ Dashed vertical lines mark the ±3<r clipping edges, 
according to our iterative procedure, as devised in Sec. 3.2. Mag residuals 
are in the sense "obs" - "syn". After outliers rejections, the global sample 
of 458 measurements has, on average, cr(Amag) = ±0.095 (see Table|7J. 



Table 7. Magnitude residuals between observed and theoretical 
magnitudes'"' 



Band 


(Obs - Syn 


All the clusters 

) 




Nout 


B 


-0.137 


0.113 


63 


10 


V 


-0.091 


0.094 


72 


1 


Rc 


+0.205 


0.065 


6 




Ic 


+0.244 


0.085 


38 


2 


J 


-0.093 


0.069 


88 


9 


H 


-0.012 


0.102 


94 


3 


K 


+0.197 


0.102 


97 




total 


0.000 <b) 


0.095 


45 S 


25 



(a) Mag residuals are in the sense of observed - synthetic one 

(b) Weigthing with the number of entries, AT; n . 

giants) No firm conclusions can be drawn, on the contrary, for the 
other seven cases, although it is evident even from a color check 



4 Note, on the other hand, the counter-example o f star #2 1 56 in NGC 679 1 , 
know n as Irr variable V70 = SBG 2240 (Moc heiska. Stanek. & Kaluznvl 
120031) and not a deviant in our spectroscopic observations. 
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Table 8. Standard synthetic photometry from SED of target stars in globular 
clusters M 71, M 15, and M 2< a ) 



Table 9. Standard synthetic photometry from SED of target stars in open 
clusters NGC 188 and NGC 6791< a ' 



M 15 



NGC 188 



ID 




V 


u 

>vc 


I 




H 




21300002+1209182 


15.1 1 


14 


.30 


13.81 


13.27 


12.58 


12 


.1)8 


1 1 .95 


21295705+1208531 


14.30 


13 


.35 


12.79 


12.18 


11.41 


10 


.82 


10.74 


21295532+1210327 


15.34 


14 


.43 


13.85 


13.24 


12.50 


1 1 


.94 


11.72 


21300090+1208571 


14.04 


12 


.90 


12.21 


11.47 


10.51 




8 1 


9.50 


21295473+1208592 


14.89 


13 


.80 


13.19 


12.55 


11.62 


1 1 


.02 


10.95 


21300461+1210327 


15.01 


13 


.85 


13.24 


12.62 


1 1.86 


1 1 


.16 


10.98 


TIIOWllMi IT IT/IT) 




13 


.46 


12.86 






10 


.79 




21300514+1210041 


15.20 


14 


.31 


13.79 


13.23 


12.45 


1 1 


,89 


11.64 


21295836+1209020 


14.80 


13 


.83 


13.27 


12.65 


11.81 


1 1 


.22 


11.01 


21295618+1210179 


14.02 


12 


.89 


12.19 


11.48 


10.49 


9 


.82 


9.55 


21295739+1209056 


14.89 


13 


.83 


13.22 


12.58 


11.65 


11 


,05 


11.03 


21300097+1210375 


14.91 


13 


.85 


13.24 


12.63 


11.85 


1 1 


.19 


11.14 


21300431+1210561 


14.76 


13 


.59 


12.91 


12.25 


11.40 


10 


,81 


10.65 


21301049+1210061 


14.45 


13 


.36 


12.73 


12.09 


11.25 


10 


.60 


10.35 


21300739+1210330 


15.22 


14 


.05 


13.36 


12.67 


11.86 


1 1 


.24 


11.10 


21300569+1210156 


15.02 


14 


.11 


13.56 


12.97 


12.21 


11 


,61 


11.50 


21300553+1208553 


15.82 


14 


.84 


14.26 


13.61 


12.50 


1 1 


.86 


11.65 


21295756+1209438 


13.80 


12 


.45 


11.72 


11.01 


10.10 


9 


.47 


9.32 


21295082+1211301 


14.61 


13 


.52 


12.91 


12.27 


11.40 


10 


.82 


10.66 


21295881+1209285< 6 > 


14.62 


13 


.40 


12.70 


12.03 


11.17 


10 


.57 


10.42 


21295716+1209175 


13.96 


13 


.03 


12.45 


11.84 


10.98 


10 


.36 


10.20 








M2 
















V 


R,, 


Id 




H 




21333827-0054569 


13.71 


12 


.72 


12.17 


11.58 


10.55 


9 


,85 


9.73 


21333095-0052154 


14.74 


13 


.68 


13.11 


12.53 


11.59 


10 


.96 


10.87 


21332468-0044252 


15.71 


14 


.65 


14.07 


13.48 


12.59 


12 


,01 


11.94 


21331771-0047273 


14.30 


13 


.1)2 


12.32 


11.66 


10.69 


9 


.94 


9.91 


21331723-0048171 


15.06 


13 


.66 


12.91 


12.21 


11.18 


10 


.44 


10.32 


21331790-0048198 


15.87 


14 


.74 


14.17 


13.54 


12.00 


11 


,05 


10.90 


21331854-0051563 


14.88 


13 


.89 


13.32 


12.72 


11.80 


1 1 


17 


11.06 


21331948-0051034 


15.11 


14 


.12 


13.55 


12.93 


12.01 


1 1 


.38 


11.17 


21331923-0049058 


16.49 


15 


.21 


14.58 


13.91 


12.52 


1 1 


,67 


11.46 


21332588-0046004 


15.42 


14 


.44 


13.89 


13.33 


12.40 


1 1 


.78 


11.59 


21333668-0051058 


14.35 


13 


.12 


12.47 


11.81 


10.76 


10 


.05 


9.92 


21333520-0046089 


14.58 


13 


.36 


12.70 


12.04 


11.02 


10 


.34 


10.23 


21333488-0047572 


14.98 


13 


.67 


12.98 


12.30 


11.33 


10 


,64 


10.42 


21333593-0049224 


15.31 


13 


.93 


13.18 


12.48 


11.49 


10 


,82 


10.60 


21333432-0051285 


14.89 


13 


.72 


13.09 


12.45 


11.50 


10 


.87 


10.77 


21332531-0052511 


15.26 


14 


.11 


13.51 


12.90 


11.98 


1 1 


.38 


11.18 


21333109-0054522 


15.98 


14 


.57 


13.80 


13.08 


12.15 


1 1 


,56 


11.37 


21333507-0051097 


15.94 


14 


.78 


14.14 


13.52 


12.63 


12 


.07 


12.02 










M71 












ID 


B 


V 


R,. 


t c 


J 


H 


K 


19535325+1846471 


14.02 


12 


.37 


11.42 


10.43 


8.92 


8 


.36 


7.93 


19534750+1846169 


14.48 


13 


.11 


12.36 


11.59 


10.41 


9 


.71 


9.57 


19535150+1848059 


13.99 


12 


.32 


11.34 


10.38 


9.10 


8 


.30 


8.02 


19535064+1849075 


14.46 


13 


.02 


12.22 


11.42 


10.27 


9 


.49 


9.25 


19534575+1847547 


14.21 


12 


.48 


11.45 


10.43 


9.12 


8 


.30 


8.02 


19534827+1848021 


14.04 


12 


.36 


11.40 


10.47 


9.24 


8 


.41 


8.09 


19534656+1847441 


14.85 


13 


.63 


12.93 


12.21 


11.22 


10 


.56 


10.47 


19535369+1846039 


15.50 


14 


.54 


14.01 


13.40 


12.56 


12 


,02 


11.85 


19534905+1846003 


14.73 


13 


.45 


12.75 


12.01 


10.91 


10 


.26 


10.15 


19534916+1846512 


14.87 


13 


.65 


13.01 


12.30 


11.19 


10 


.46 


10.30 


19534178+1848384 


15.77 


14 


.55 


13.87 


13.17 


12.24 


1 1 


.58 


11.40 


19535676+1845399 


15.66 


14 


.50 


13.86 


13.18 


12.22 


1 1 


.58 


11.43 


19533962+1848569 


15.74 


14 


.70 


14.13 


13.50 


12.50 


1 1 


.87 


11.74 


19533864+1847554 


15.39 


14 


.1(1 


13.42 


12.72 


11.72 


11 


.19 


11.09 


19534615+1847261 


14.68 


13 


.21 


12.38 


11.55 


10.35 


9 


.60 


9.44 


19535610+1847167 (c) 


14.96 


13 


.27 


11.21 


9.26 


7.89 


7 


.08 


6.83 


19534941+1844269( d > 


13.88 


12 


.1)2 


10.71 


9.42 


7.96 


7 


.20 


6.96 



(a) After correction for the systematic offsets, according to Table|7] 

(b) Dropped: I outlier 

(c) Dropped: SR vari ahle 7 Sge- R V outlier; 

(d) Var AN 48.1928 i Baa de 1928 



ID 


13 


y 


D 

K C 


T 

L C 


j 


^1 


1^ 


00445253+8514055 


13.65 


12.34 


11.66 


11.01 


10.04 


9.41 


9.45 


00475922+8511322 


13.68 


12.13 


11.38 


10.73 


9.78 


9.20 


9.25 


00465966+8513157 


13.87 


12.40 


11.71 


11.11 


10.25 


9.7] 


9.67 


00453697+85 15084< 6 ' 


15.49 


13.90 


13.10 


12.40 


11.44 


10.87 


10.77 


00442946+8515093 (c) 


15.84 


14.28 


13.44 


12.73 


1 1.81 


11.27 


11.32 


00473222+851 1024 (d) 


15.79 


14.24 


13.40 


12.73 


11.87 


11.40 


11.49 


00554526+8512209 


12.25 


10.82 


10.06 


9.34 


8.31 


7.59 


7.54 


00463920+8523336 


12.89 


11.56 


10.87 


10.21 


9.19 


8.53 


8.50 


00472975+8524140 


14.05 


12.91 


12.33 


11.74 


10.70 


10.11 


9.98 


00441241+8509312 


12.81 


11.28 


10.45 


9.68 


8.60 


7.85 


7.78 


00432696+8509175 


14.17 


13.14 


12.61 


12.08 


11.25 


10.76 


10.59 


0047 1847+85 19456 (e) 


14.89 


13.03 


12.19 


11.47 


10.54 


9.99 


10.07 


00461981+8520086^ 


14.58 


12.50 


11.60 


10.88 


9.98 


9.43 


9.40 


00463004+8511518 (b) 


15.67 


14.20 


13.45 


12.78 


11.91 


11.42 


11.38 


00490560+8526077 


13.89 


12.65 


12.02 


11.41 


10.48 


9.90 


9.97 


00420323+8520492 


11.29 


9.75 


8.91 


8.14 


7.05 


6.43 


6.20 








NGC 6791 










ID 


B 


V 


R,. 


Io 


J 


H 


K 


19210807+3747494 


15.29 


13.98 


13.33 


12.70 


11.66 


11.06 


10.87 


19204971+3743426 (fe) 


15.52 


13.97 


12.31 


10.61 


9.04 


8.20 


7.97 


19205259+3744281 


15.69 


14.10 


13.19 


12.32 


11.17 


10.45 


10.20 


19205580+3742307 


16.23 


14.90 


14.23 


13.59 


12.61 


12.03 


11.96 


19205671+3743074 


15.97 


14.64 


13.96 


13.31 


12.32 


11.74 


11.66 


19210112+3742134 


15.94 


14.48 


13.68 


12.92 


11.82 


11.06 


10.96 


19211606+3746462 


15.29 


13.74 


12.11 


10.45 


8.9] 


8.10 


7.84 


19213656+3740376 


15.58 


14.07 


13.24 


12.46 


11.40 


10.72 


10.50 


19210326+3741190 


15.72 


14.39 


13.69 


13.01 


12.03 


11.44 


11.36 


19213635+3739445 


16.10 


14.76 


14.1 1 


13.49 


12.47 


11.78 


11.58 


19212437+3735402 


15.83 


14.46 


13.79 


13.15 


12.14 


11.45 


11.30 


19212674+3735186 


15.33 


13.98 


13.25 


12.56 


11.59 


10.99 


10.79 


19211632+3752154 


15.22 


13.98 


13.32 


12.67 


11.71 


11.12 


10.93 


19211176+3752459 


15.65 


14.38 


13.73 


13.11 


12.13 


11.47 


11.41 


19202345+3754578 


14.57 


12.80 


11.43 


9.63 


7.96 


7.08 


6.84 


19205149+3739334 


13.26 


11.66 


10.23 


8.75 


7.35 


6.58 


6.23 


19203285+3753488 


14.97 


13.39 


1 1.67 


10.00 


8.46 


7.55 


7.19 


19200641+3744452 


14.67 


13.34 


12.61 


11.89 


10.79 


10.06 


9.87 


19200882+3744317 


15.17 


13.62 


11.61 


9.67 


7.92 


7.07 


6.72 


19203219+3744208<*> 


16.22 


14.77 


12.46 


10.23 


8.18 


7.29 


6.93 



(a) After correction for the systematic offsets, according to Table|7] 

(b) Dropped: B, J outlier; 

(c) Dropped: B, J, H outlier; 

(d) Dropped: B, J outlier; 

(e) Dropped: B, J outlier; 

(f) Dropped: B, J, H outlier; 

(g) Dropped: B, J outlier; 

(h) Dropped: I outlier; V13 - Var7 Ide. Marchi et al l?.007n 

(i) V70 = SBG 2240: Irr var (Mocheiska, Stanek, & Kaluzny 2003) 



(see Fig.|4]l that they have been picked up in an intrinsically differ- 
ent status with respect to previous data in the literature0 We there- 
fore commend the stars in the list of TablefTOlas special candidates 
for further in-depth investigations for variability. 



4 SPECTRAL ENERGY DISTRIBUTION AND 
BOLOMETRIC LUMINOSITY 

The synthetic photometric catalogs obtained from the observed 
spectral database had a twofold aim: firstly, this procedure allowed 



5 Curiously enough, however, one may note that 6 out of the 7 remaining 
objects are all located in NGC 188, and are both B and J outliers. Both 
photometric bands actually cover the "bluer" wavelength regions of both 
the LRS and NICS spectra, respectively. This coincidence might perhaps 
indicate some hidden problem with the flux calibration procedure during 
the observation of this cluster. 
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Figure 4. Color distribution of photometric outliers, according to our 3cr 
clipping procedure (see Fig.[5J- Target location for the whole star sample in 
the synthetic vs. observed color planes are displayed, with dark solid dots 
marking the "dropped" objects (see Table llOt . 



Table 10. Stellar outliers of our sample in 


the different photometric bands 


Cluster 


ID 




Outlier in 


Notes 






B V 


Rc Ic J 


H K 


M15 


ISO 




X 




M7I 


4212 


X X 




ZSge-SRaP~ 175 d 




5755 






V2 = AN 48.1928 - Ir 


N188 


567 


X 


X 






(,52 


X 


X 






1109 


X 


X 






630 


X 


X 






174 


X 


X 






1352 


X 


X 




N6791 


35(12 




X 


V13-Var? 


Tot 




7 1 


3 6 


2 



us to self-consistently match broad-band magnitudes with the in- 
ferred measure of rriboU m order to obtain the corresponding value 
of the bolometric correction; secondly, the study of the magnitude 
residuals with respect to the literature data provided us with the ap- 
propriate offsets in flux rescaling such as to "smoothly" connect our 
optical and infrared spectra and lead therefore to a more accurate 
estimate of mboi- 

Operationally, for the latter task, we proceeded as follows. 
Taking into account the individual set of (Obs — Syn) magnitude 
residuals, for each star in our sample we computed a mean optical 
and infrared offset (Aitilrs and ArriNics, respectively) by sepa- 
rately averaging the B ,V, R c , I c and J, H, K mag residuals. The 
LRS spectra and the NICS observations have then been matched by 
multiplying visual and IR fluxes by a factor 10 -0 ' 4 ( Am LRs) anc j 
^q-0.4 (Am NICS )^ reS p ec ti ve iy Foreground reddening has been 
corrected relying on the standard relation fc(A) = A{X) /E(B — V) 
(Scheffler 2006), where the appropriate value of the color excess 
E(B — V) is from the headers of Tables Q] to [5] In its final form, 
the SED is reshaped such as / D (A) = /(A) 10 4 k WB(B-V)_ 
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Figure 5. The resulting (dereddened) SED according to optical and infrared 
observations for an illustrative stellar subset of each cluster, including the 
brightest (and roughly coolest) and faintest (i.e. warmest) stars. Note, espe- 
cially for the M 15 stars, the strong impact of telluric water vapor bands at 
1.38 and 1.88 fj,m. Their variability along the observing nights prevented, 
in some cases, any accurate cleaning procedure. See discussion in Sec. 3.1. 
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The LRS and NICS spectra have been connected at 8800 A, 
by smoothing the wavelength region between 7800 and 10000 A 
(in order to gain S/N, especially for LRS poor signal at the long 
wavelenegth edge). In Fig. [5] we summarize our results for an illus- 
trative set of SEDs by picking up for each cluster the brightest (i.e. 
roughly the coolest) and faintest (i.e. warmest) stars in our sampleQ 
Note, fr om the figure, the s t ricking presence of the CO bump about 
1.6 fim dFrogel et al ] |l978l : lLancon & Mouhcine 2002), as well as 
the broad H2O absorption bands to which the sharper (and variable) 
emission of telluric water vapor superposes (see, in particular, the 
case of M15 stars in the figure). This made far more difficult any 
accurate cleaning procedure, as we discussed in Sec. 3.1. 



4.1 Temperature scale 

Although sampled over a wide wavelength range, SED of our stars 
still lacks the contribution of ultraviolet and far-infrared luminosity. 
Clearly, a safe assessment of this contribution is mandatory to lead 
to a confident measure of the bolometric magnitude. As the amount 
of energy released outside the spectral window of our observations 
critically depends on stellar temperature, our task to compute BC 
requires in fact a parallel calibration of T c s in the range of our red 
giant stars. 

Among the many outstanding efforts in this direction. 



and 



1979) 



1980) 



( 2000), VandenBerg & Clem 
Worthev&Leei d2006h. 



we have to recall the works of Flower! dl975l). iBesselll 
Blackwell, Petford. & Shalliiil dl980h. iRidgwav et all 
Bessell. Castelli. & Pled dl998l>. iHoudashelt. Bell. & Sweigartl 

( 120031) . iBertone et al.l ( bOM) 
In their exhaustive analysis, 
Al onso, Arribas, & Martinez-Roger ( 1999) provided an ac- 
curate analytical set of fitting functions, that calibrate stellar 
effective temperature vs. Johnson/C o usins broad-band colors. The 
lAlonso. Arribas. & Martmez-Roged dl999h calibration relies on 
the IRFM estimate of stellar surface brightness, and considers stars 
of spectral type K5 or earlier, spanning a wide metallicity range 
(-3.0 < [Fe/H] < +0.2). Within this range, the Authors claim 
an internal accuracy in the definition of T a g better than 5%. As a 
further important result of their work, some colors, like (V — /), 
(V — L'), (J — K) and (/ — K) are found to be fair tracers of 
temperat ure, almost independently from stellar metall icity. 

The lAlonso, Arribas. & Martinez-Roger (1999) calibration, 
however, strictly applies only to stars warmer than ~ 4000 K, 
while our stellar sample definitely spans a wider color range. This 
is certainly the case, for instance, of the brightest giant stars in 
NGC 6791, too (infra)red to match the Alonso et al. fitting func- 
tions. For these cases one could rely on the wider validity range 
of the (B — V) calibration, although the advantage may only be a 
nominal one as any optical color, like (B — V) tends na turally to 
saturate when moving to r e g < 4000 K dJohnsonlll966l see also 
Fig. 2 in lAlonso. Arribas. & Martmez-Rogej|l999l) . 

Considering the whole set of the Alonso et al. fitting functions, 
we eventually chose four reference colors to assess the value of 
effective temperature for our stars. Two colors, namely (B — V) and 
(J — K) are entirely comprised within the LRS and NICS spectral 
branches, respectively, and they can therefore ostensibly probe the 
shape of SED in a more self-consistent way. To these two colors 
we also added (V — I c ) and (V — K), as they provided a check of 



6 For the interested reader, the entire spectral database is available in elec- 
tronic form upon request, or directly on the web at the authors' web site 

|http : / /www .bo .astro . it /^eps/home .html| 



60 



40 - 



20 - 



PL 



332 meas. 
a(T)=150 K 



-1000 



-500 500 

T-<T> 



1000 



Figure 6. Histogram of temperature difference for all the data reported in 
Tables[8]and|9](cols. 6 to 9) with respect to the adopted mean estimate ((T) 
of col. 10). A total of 322 entries are available for the whole stellar sample. 
The resulting distribution gives a direct measure of the internal uncertainty 
of our temperature scale, amounting of o"(T ef j) = ±150 K for the standard 
individual estimate. 



our flux calibration bridging the optical and infrared regions of the 
spectra. 

Dereddened colors for each stars in our sample provided even- 
tually a set of nominal values of T e g, by entering the appropriate 
fitting functions. The "allowed" values of T B g (i.e. if comprised 
within the boundary limits of the adopted calibration functions) 
were then averaged, deriving the mean fiducial value of the effec- 
tive temperature, reported in Tablell llandll2l(column 10). In case 
of just one T e g estimate (typically from (J — K) color) we also 
added the (V — K) output (repor ted in italics in the tables) trusting 
on a fa irly smooth trend of the lAlonso. Arribas. & Martmez-Rogen 
calibration for this color, when extrapolated to cooler tem- 
peratures (see Fig. 8 and Fig. 10 therein). 

Once combining the different temperature estimates from the 
four reference colors in our analysis, we report in Fig.|6]the result- 
ing T — (T) distribution, considering the whole set of 322 individ- 
ual residuals. The figure confirms that an unbiased estimate of T c g 
may eventually be achieved with our procedure, within a ±150 K 
uncertainty on the standard measure. As, typically 2-4 useful tem- 
perature estimates are available from the colors of each star (see, 
again, Table[TT1andll2t. we may expect final T c s values for our 
sample to be assessed within a 70-100 K (i.e. 1-3%) internal uncer- 
tainty. 



4.2 Toward m bo i 

The fiducial effective temperature, as reported in col. 10 of TablefTTI 
and[l2] provided the reference quantity to constrain the unsampled 
fraction of stellar luminosity, outside the wavelength limits of our 
spectral observations. No univocal procedure can be devised to ef- 
fectively tackle this problem; from one hand, in fact, both the ul- 
traviolet and mid- and far-infrared stellar emission can in principle 
be modulated by a number of different mechanisms (mass loss and 
stellar winds, or circumstellar gas and dust lanes thermalizing ultra- 
violet and optical photons, photospheric spots, pulsating variability 
etc.). On the other hand, one would better like to proceed with a 
straight heuristic approach, such as to self-consistently size up the 
amount of "overflown" luminosity and decide the accuracy level in 
its correction procedure, according to an "ex-post" analysis of the 
results. 
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Table 11. Inferred temperatures, bolometric magnitude and bolometric corrections for target stars in globular clusters M 71, M 15, and M 2 
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M 15 



ID 


(B-V)„ 


(V-I c )„ 


(V-K) D 


(J-K) 


T BV 
"K 


T v] 
°K 


T VK 
"K 


T JK 

"K 


CO 

°K 


Bol 


BC V 


BC A - 


21300002+1209182 


0.71 


0.87 


2.08 


0.58 


5043 


5030 


5010 


4748 


4962 


13.742 


-0.25 


1.83 


21295705+1208531 


0.85 


1.01 


2.34 
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4399 
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4660 
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3941 


4213 


11.934 


-0.65 


2.42 


21295739+1209056 


0.96 


1 .09 


2.53 


0.57 


4574 


4549 


4549 


4782 


4614 


13.068 


-0.45 


2.07 


21300097+1210375 


0.96 


1.06 


2.44 


0.66 


4574 
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4618 


4490 
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-0.41 
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0.81 
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0.66 
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4694 
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-0.33 
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21300553+1208553 
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1.07 


2.92 
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4680 


4588 


4249 


4136 


4413 


14.012 


-0.52 


2.40 


21295756+1209438 


1.25 


1.28 


2.86 


0.73 


4230 


4227 


4291 


4307 


4264 


11.556 


-0.58 


2.27 


21295082+1211301 


0.90 


1.09 


2.59 


0.60 


4536 


4540 


4498 


4414 


4499 


12.778 


-0.43 


2.15 
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21333827-0054569 


0.93 


1.04 


2.83 


0.70 


4640 


4639 


4312 
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4437 


11.990 


-0.54 
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9 1 33300 s , 00^9 1 ^4 


I 00 


1 ()9 


9 


n 60 


4940 


4 6 1 8 


44 47 


4412 


4504 


1 1 ()99 


—0 47 


9 1 8 


21332468-0044252 


1.00 


1.07 


2.55 


0.62 


4540 


4570 
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4566 
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-0.43 
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21331771-0047273 
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1.26 


2.95 


0.75 


4251 


4250 


4230 
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4246 
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2.30 
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1.35 


3.18 
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12.721 
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2.42 
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3.68 
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3912 
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4017 
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0.93 


1.07 


2.67 


0.71 


4640 


4579 


4431 


4358 
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13.231 


-0.47 


2.10 


21331948-0051034 


0.93 


1 .09 


2.79 


0.81 


4640 
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4340 


4111 
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13.428 
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2.28 
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14.083 


-0.94 
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1.01 
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4655 


4701 


4416 
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2.25 
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1.21 
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0.81 
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2.37 
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1.22 


2.97 


0.76 
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4312 


4217 


4229 
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12.528 


-0.65 


2.32 
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1.25 


1.27 


3.09 


0.88 


4215 
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4143 


3959 
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12.813 


-0.67 


2.41 
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1.35 


3.17 
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2.42 
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0.70 
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19535325+1846471 


1.39 


1.54 


3.75 


0.86 
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3881 
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2.61 
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11.760 
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0.96 


1.02 


2.47 
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0.58 
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0.96 


0.95 


2.66 


0.76 


4710 


4849 


4447 
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12.980 
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1.21 


1.26 


3.08 
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On this line, we therefore decided to proceed in the most 
straightforward way for each star, by extrapolating its observed 
SED to both ultraviolet and infrared windows by means of two 
black-body branches, of appropriate (fixed) temperature (T) as 
in Tables [TT] and [12] The two spectral branches have been sep- 
arately rescaled to the (dereddened) flux values of the observed 
SED by setting the boundary wavelengths respectively at 4000 and 
22500 A; the integrated luminosity has then been computed within 



the three relevant regions of each stellar SED, identifying the ultra- 
violet contribution Ztjv (between ^ A ^ 4000 A) an optical/mid- 
infrared luminosity Z b s (4000 ^ A ^ 22500 A) and a far-infrared 
contribution Zpir (longward of 2.25 /jm). For comparison, the same 
excercise has been repeated for a straight black-body spectral dis- 
tribution exploring the luminosity fraction emitted shortward of 
A ^ 4000 A and longward of A > 22500 A along the tempera- 
ture range of our sample. 
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Table 12. Inferred temperatures, bolometric magnitude and bolometric corrections for target stars in open clusters NGC 188 and NGC 6791 



NGC 188 
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4478 


5078 


4460 


11.311 


-0.56 


2.09 


00465966+851315 
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4078 
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0.93 
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0.62 
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2.00 


00420323+852049 


1.45 


1.48 


3.32 
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Our results are summarized in Fig. [7] Compared to the black- 
body aproximation, real stars are brighter at longer wavelength and 
slightly fainter, on the contrary, at UV wavelength. In total, one 
sees from Fig.[7jthat the fraction of "lost" luminosity, namely F\ — 
(Iboi — fcbsV^boi, turns to be about 15% for the bulk of red giants in 
our sample; this figure can however quickly raise with decreasing 
temperature, and about 1/3 of bolometric luminosity might in fact 
be "stored" at FIR wavelengths. Within these limits, and accounting 
for the 70-100 K internal uncertainty of our temperature scale, one 
sees from the Fig. [JJthat ra^ \ can be secured for our sample stars 
within a few 0.01 mag uncertainty^ 

Starting from the bolometric flux (which also includes the un- 
sampled luminosity fraction, according to our procedure), the ap- 
parent magnitude for each star derives as iriboi = — 2.5 log /boi + 
Z.P.. If we assume for the Sun an absolute Mboi = +4.72, and 
Lq = 3.89 10 33 ergs -1 , the bolometric zero point directly de- 
rives as Z.P. = —11.50 mag. On the same line, the BC scale is 
fixed once adopting an obs erved value f or the apparent V magni- 
tude of the Sun. Following Hangl Jl99lh . if m® = -26.78, then 
M§ = +4.79 and a BC$ = -0.07 mag derives. Our output, for 
the whole stellar sample, is reported in col. 1 1 of Tablell 1 landl 121 
together with the relevant (dereddened) BC to the V and K photo- 



7 The claimed m^ i uncertainty simply derives as a ~ 8Fi/dT e f{ X 
<x(T e fj), where cr(T e ff) < 100 K and the F\ derivative can be estimated 
from Fig. [7j In any case, it is clear from the figure that, by neglecting any 
further luminosity correction to our data for the unsampled luminosity, we 
would be overestimating m^oi by at most 0.3 mag. 




5000 4000 3000 
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Figure 7. Estimated fraction of unsampled stellar luminosity for the stars in 
our sample (big solid dots). The relative contribution to stellar bolometric 
luminosity from lost emission at short (i.e. for A 4000 A, small square 
markers on the plot) and long (i.e. for A > 2.25/^m, small triangles) wave- 
length is sized up by extrapolating the observed SED with two black-body 
(BB) "wings" at fixed (T), as from col. 10 of Table [Til and 1 121 The same 
excercise is carried out for a full BB spectrum along the 5500-3000 K tem- 
perature range (dashed lines labelled "UV" and "IR" for the short and long 
wavelength contribution, respectively, together with their summed contri- 
bution, as in the solid line). Compared to a plain BB case, note that real 
stars at cool temperatures display a brighter IR luminosity. 
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Figure 8. The BC vs. color (left panels) and BC vs. T c ff (right panels) distribution of our stellar sample (dots and triangles, for metal-poor and metal-rich 
stars, respectively). Synthetic colors have been corrected for Galactic reddening. Solid lines are our derived calibrations, according to the set of eqs. (4) and 
®. 



metric bands (BCv and BCk, respectively in col. 12 and 13 of the 
tables). 



5 RESULTS AND DISCUSSION 

The data of Tables [TT1 and [121 are the main output of our analysis. 
According to our results, we can explore three relevant relation- 
ships, linking BC with the effective temperature of stars and with 
two reference colors like ( B — V ) and (V — if). Given the tempera- 
ture range of red giants, it could be of special relevance to consider 
the iv-band BC; however, for its more general interest, we will also 
include in our discussion the more standard case of the BCv. 



5.1 BC-color-temperature relations 

Like for a color-color diagram, the BC vs. color relationship can 
be regarded as an intrinsic (i.e. distance-independent) feature char- 
acterizing the stellar SED. On the corresponding theoretical side, 
we want also to study here the resulting dependence of BC on stel- 
lar effective temperature, a relation that allows us to more directly 
match the observations with the theoretical predictions of stellar 
model atmospheres. 

In a first set of plots (see Fig. [8), we display the observed 
distribution of our stars in the different planes. In order to single 
out any possible dependence on chemical composition of stars, we 



marked differently metal-poor ([Fe/H] < —1.0 dex, dots) and 
metal-rich ([Fe/H] > — 1.0 dex, triangles) objects. For better con- 
venience in our study, we also fitted the overall distribution analyt- 
ically; a useful set of fitting functions for the BC vs. T e fr relations 
along the 3300 < T ff < 5000 K temperature range results: 



BC V = -exp(27500/T off )/1000 
BC K = -6.751og(T cff /9500) 



(0.11,0.989) 



= (0.05,0.978). 



(4) 



As for the color relations, the non-monotonic trend of BCv 
vs. (B — V) (see left upper panel in Fig. [8) prevents us to use the 
color as independent (i.e. "input") variable in our fit. In this case we 
had therefore to adjust an inverse relation, assuming BC as the run- 
ning variable. The corresponding set of analytical solutions, along 
the same temperature range of the previous equation set, eventually 
results: 



B - V = 1.906 BC v exp(BC\ 



(obv.p) =(0.11,0.863) 



(5) 



V-K = 1/(1 - 0.283 BC K ) 

(tvk,p) =(0.13,0.991), 

All these fits are superposed to the data of Fig.[8]as a solid line, 
lust on the basis of our data note how difficult it is to firmly 
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Figure 9. The distribution of BC residuals for our stellar sample vs. cluster 
metallicity. The displayed ABC is intended as the difference between the 
values of cols. 12 and 13 of Tables [TTI and [T2l and the output of eq. @ 
entering with the adopted effective temperature of stars, as in col. 10 of the 
tables. Note the lack of any evident correlation with [Fe/H], as discussed in 
more detail in Sec. 5.2. Data in the plot have been slightly spread around the 
cluster [Fe/H] value for better reading. Dot size is inversely proportional to 
star temperature (i.e. bigger dots = cooler red giants). 



constrain the (B — V) vs. BCv behaviour at very low temperature. 
From one hand, in fact , the intervening effec t of the TiO absoprtion 
at visual wavelength dKucinskas et al.l 2005) makes the (B — V) 
color of stars cooler than ~ 3700 K to strongly saturate reaching a 
maximum of about (B — V) max — 1.5 and turning back to bluer 
values for later M-type stars. On the other hand, the apparent trend 
of our sample in this range is evidently biased by the NGC 6791 
stellar population with just a few super metal rich giants constrain- 
ing the BCv trend at the most extreme negative values. 



5.2 BC response to metallicity 

As a part of our observing strategy, the sampled stellar population 
of the five clusters would in principle allow to better single out 
any possible dependence of BC on stellar chemical composition. 
As far as Helium content is concerned, for insta nce, this problem 
has already been tackled by Girardi et al. (2007) through a series 
of theoretical models based on the lKuruczl l l 19921) ATLAS9 model 
atmospheres. As a main result of their discussion, these authors did 
not find any relevant impact on stellar BC to optical photometric 
bands when Helium changes up to AY = +0.2, for fixed effective 
temperature. To some extent, this is a not so surprising behaviour; 
Helium is in fact a substantial contributor to mean particle weight 
of stellar plasma but a negligible contributor to chemical opacity. 
Accordingly, with varying Y in the chemical mix, one has to expect 
a much more explicit impact on stellar temperature for fixed mass 
of stars, rather t han on colors or SE P for fixed effective temperature 
(as explored bv lGirardi etai]|2007l models, indeed). 
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Figure 10. Temperature distribution of red-giant stars in each of the five 
clusters of our sample, according to Tablell 1 land l 121 Line tickness is pro- 
portional to the star density along the spanned temperature range. Note that 
only cluster NGC 679 1 contains stars cooler than ~ 3800 K. 



The situation might in principle be different for the metals, 
mainly through their pervasive effect on stellar blanketing at short 
wavelength. In addition, metals are the basic ingredients required 
to produce molecules like TiO, SiH or CH, whose impact may 
be extremely relevant at blue and visual waveleng th, when ef- 
fective temperature lowers below 3500 K ( Kucinskas et alj|2005l ; 
iBertone et ai1l2008h . 

Taking the results of Tablesl 1 1 land[L2las a reference, in Fig. [9] 
we plot the BC residual distribution computed as a difference be- 
tween the inferred BC (cols. 12 and 13 in the tables) and the "mean" 
locus of eq. (0), once entering the equations with the fiducial (T) of 
col. 10. The BC residuals are displayed along the [Fe/H] distribu- 
tion of the five star clusters, as labelled on the plots. Just a glance to 
both panels of the figure makes evident the lack of any drift of BC 
with stellar metallicity. Within the accuracy limits of our analysis, 
this means that two red giant stars of the same effective temperature 
but different [Fe/H] have virtually indistiguishable values of BC to 
V and K bands. 

On the other hand, to correctly understand our conclusion, 
one has to pay attention to the different temperature regimes that 
mark spectral properties of red-giant stars. In fact, stars warmer 
than ~ 4000 K may have their SED depressed at short wavelength 
mostly in force of atomic transitions of Fe and other metals; on 
the contrary, for a cooler temperature, the metal opacity mainly 
acts in the form of molecular absorptions, making the broad band 
systems the prevailing features that modulate the stellar SED. As 
a consequence, while for stars of spectral type G or earlier any 
change of Z simply implies a change in the blanketing strength, 
this may not straightforwardly be the case for later spectral types, 
where molecules play a much more entangled role with changing 
Tiff- 
in order to better quantify the terms of our analysis, in this 
respect, we display in Fig.[l0]the temperature distribution of stars in 
our sample across the metallicity range spanned by the five clusters 
considered. As a stricking feature, note that only for NGC 679 1 we 
are able to probe stars cooler than ~ 3800 K. The obvious caveat 
in our discussion is therefore that we can only assess the impact of 
atomic blanketing on stellar BC, while no firm conclusions can be 
drawn for the BC dependence on molecular absorption, facing the 
evident bias of our star sample against cool (T e a <C 4000) objects. 

As far as the blanketing is the prevailing mechanism at work 
in G-K stars, basic physics of stellar atmospheres leads to con- 
clude that the V-band (and even more the K-bwiA) luminosity are 
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Figure 11. Residual (B — V) and (V — K) distribution vs. stellar tempera- 
ture for stars in Table ll 1 l and l 121 Color residual is computed as a difference 
between observed and expected values by entering eq. (3} with the BC out- 
put of eq. @. Metal-poor and metal-rich stars are singled out by diamonds 
and dot markers, respectively, taking the value [Fe/H] = —1.0 dex as a 
reference threshold. Marker size increases with [Fe/H], throughout. 



nearly unaffected by metal absorption, so that BC cannot vary much 
with [Fe/H]. Rather, B (and even more U) magnitudes must be 
more strongly modulated by metal abundance making BCs (and 
BCc/) more directly sensitive to [Fe/H]. On the other hand, as 
BCb = BCv — (B — V), one can straight "translate" this metal- 
licity effect in terms of apparent (B — V) color change. This is 
shown in Fig. QT| where for each star in our sample we computed 
the residual (B — V) and (V — K ) color as a difference between ob- 
served and expected values by entering eq. ((5) with the fitted value 
of BC as from eq. {4j. Metallicity is traced in the plot by the marker 
size (the bigger the marker the higher the [Fe/H] value); again, we 
discriminate between metal-poor (diamonds) and metal-rich (dots) 
stars, taking the value [Fe/H] = —1.0 dex as a reference thresh- 
old. 

A trend of A(B — V) vs. cluster metallicity is now clearly ev- 
ident, with the metal-poor and metal-rich star samples neatly seg- 
regated in the plot, the latter stars displaying a "redder" (B — V) 
color (and correspondingly a positive color residual) for fixed ef- 
fective temperature. On the contrary, note that both "metal-poor" 
and "metal-rich" stars are well mixed in the A(V — K) plot, wit- 
nessing once more the property of the V — K color as a virtually 
metal-independent feature. 

Considering in more detail the A(B — V) distribution vs. clus- 
ter metallicity, a fit to the data provides 



-ABCb = A(B~V) = 



0.10 

±1 



[Fe/H] +0.13 



±2 



(6) 



with error bars at la level and (rms, p) = (0.09 mag, 0.70). 



8 Of course, following our previous arguments, we had to exclude from our 
analysis cluster NGC 679 1 , for its obvious bias in containing the emipirical 
T e ff vs. (B — V) relationship for stars at supersolar metallicity. 



5.3 Comparison with other BC scales 

For a better understanding of our results it is relevant to com- 
pare our output with other popular calibration scales often taken 
as a reference in the current literature and especially attempting 
to extend their analysis to cool (T c g < 3500 K) stellar tem- 
peratures. In particular, we will focus here on different theoret- 
ical BC calibrations relying on the three leading codes for ad- 
vanced computation of stellar model atmospheres, namely AT- 
LAS 9 dKuruczl 1 1994 hereafter la belled as "AT9"), NEXTGEN 
dHauschildt. Allard. & Baronl [19991 . "NG"), both as reported by 
iBertone et ail d2004l) . and MARCS dBell & Gustafssorj 1 19781 as 
adopted by Houdashelt et al. 20 00, "HOP' label) also in its up- 
dat ed versions (NMARCS, as in Plez. Brett, & Nordlundl d 1992b 
and lBessell. Castelli. & Plezldl998L "NM"). 

We will also co nsider in our a nalysis two empirical stud- 
ies, i.e. the ones of IJohnsonl dl96d referred to as "J66") and 
iMontegriffo et alj d 19981 labelled as "M98"), both based on a care- 
ful analysis of infrared colors to assess the problem of the bolomet- 
ric correction and a self-consistent temperature scale for red giant 
stars. All the bolometric scales in the figure have been shifted such 
as to agree with our assumption that BCy = —0.07 mag. 

A synoptic look of the different theoretical and empirical 
frameworks is eased by the four panels of Fig. [12] where we re- 
port the BCv and BCa- scales vs. observables (i.e. (B — V) and 
(V — K) colors, respectively) and theoretical (T e ff) reference quan- 
tities. In all respects, this figure is fully equivalent with, and can be 
compared to, Fig. [8] where we reported our own results. 

Just a quick look to the different curves of Fig. [12] gives an 
immediate picture of the inherent uncertainties in predicted BC ac- 
cording to the different calibration scales. The big issue, in this re- 
gard, much deals with the way models can reproduce cool stars and 
observations can account for the (B — V) "saturation" vs. temper- 
ature consequent to the shifted emission toward longer wavebands 
when stars become cooler than 3500 K. This effect makes the B- 
luminosity contribution to drop to nominal values among red gi- 
ants, and the increasingly important role of molecular absorption 
strongly modulates optical colors of K- and M-type stars. 

The still inadequate theoretical performance in modelling such 
cool stars with convenient accuracy fatally frustrates also any em- 
pirical effort to derive a firm temperature scale and an accurate 
abundance analysis for stars at the extreme edge of the temperature 
distribution (see, e.g.. IBertone et alj|2008l and lolling et al]|2009l 
for useful considerations on this subject). 

As far as the BCv vs. (B — V) behaviour is concerned, the ref- 
erence calibrations display the largest spread, with M98 predicting 
increasingly redder stars with decreasing temperature. At the op- 
posite, NM predicts a sharp color "turnback", with BCv increasing 
in absolute value among cool stars getting bluer and bluer. Def- 
initely, the emipirical calibration by J66 still remains a reference 
one, fairly well tracking the observations. This trend is very closely 
replied also by the MARCS models by H00, that provide an even 
better match to the data and a substantial agreement with our fitting 
function as in Fig. [8] 

By converting colors to the theoretical plane of effective tem- 
perature (right upper panel of Fig. l 12b . the picture slightly changes, 
in particular with a stricking discrepancy of the J66 and the theoret- 
ical NG temperature scale for T c g < 3800 K. Both sources predict, 
in fact, much shallower corrections for cool stars than we observe. 
An overall agreement has to be reported, on the contrary, among 
the other calibrations, all replying our eq. {4). 

The situation is much eased in the infrared domain, where a 
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Figure 12. Same as Fig. [8] but com paring our data with di fferent theoretical and empirical cal ibrations from Ijohnsonl j 19661 "J66" labels), Be rtone et all 
120041) using ATLAS9 jk urucz 199 ^, "AT9") and NextGen Hauschildt . Allard. & Barod fl999l "NG") synthesis codes for model atmosph ere computation, 
Montegriffo et all Jl998t "M98'").lHoudashelt. Bell & Sweigartl feOOOl ' 'HOP") using MARC S theoretical code by iBell & Gustafsson 1978), and its updated 
version (NMARCS), as in lPlez. Brett, & Nordlund Jl992h and lBessell. Castelli. & Pled J 19981 "NM"). 



monotonic relationship between (V — K) color and BCa' charac- 
terizes red giants stars. In this new framework both the theoretical 
and empirical planes are well reproduced by the different calibra- 
tion scales, with the only remarkable exception of J66 that, to some 
extent "allows" stars to store a bigger fraction of their bolometric 
luminosity in the infrared. This leads to a tipping BCk — 2.7 and 
a too "red" (V — K) for a given value of T c ff . 

Combining the different pieces of information coming from 
these comparisons, it seems that the H00 MARCS models are by 
far the best ones in matching our BC estimates, closely replying in 
every panel of Fig. [12] our empirical fitting functions of eqs. {4j, 
(|5) and Fig. [8] In spite of this comforting appearence, however, 
this conclusion may be even more puzzling from a physical point 
of view, as the H00 models have been a fortiori tuned up such as 
to reproduce the observed colors of M stars. As described by the 
authors, this required in particular to strongly enhance the assumed 
TiO opacities well beyond the admitted physical range suggested 



by mol ecular theory and impl emented in the "standard" MARCS 
library dGustafsson etai]|2008h . 



6 SUMMARY AND CONCLUSIONS 

The firm knowledge of a fully reliable link between observations 
and stellar evolution models is a basic, crucial requirement for any 
safe use of stellar clocks and population synthesis templates in the 
study and interpretation of the integrated spectrophotometric prop- 
erties of distant galaxies. Actually, the "stellar path" to cosmology 
is strictly dependent, among others, on the accurate determination 
of the bolometric emission of stars, with varying effective temper- 
atures and chemical abundance. 

In this framework, we have tackled the central question of the 
possible BC dependence on stellar metallicity by securing spec- 
troscopic observations for a wide sample of 92 red-giant stars in 
five (3 globular + 2 open) Galactic clusters along the full metallic- 
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ity range from [Fe/H] = —2.2 up to +0.4 (see Sec. 3). Spectra 
cover the wavelength range from 3500 A to 2.5^m, collecting op- 
tical and IR observations. A delicate task for the final settlement 
of our stellar database dealt with the accurate flux calibration and 
a consistent match of the optical and near-IR sides of the spectra 
such as to reproduce, for each star, the broad-band BV R C I C J HK 
photometry available in the literature (Sec. 2 and 3). Allover, we 
are confident that stellar SED along the entire sampled wavelength 
range has been set up within a ±10% internal accuracy (see Table|7J 
and Fig. [3}. 

According to our previous arguments, however, one has also to 
carefully account for the lost contribution of ultraviolet and far-IR 
luminosity to the bolometric flux, depending on the effective tem- 
peratu re of stars. Based on the lAlonso. Arribas. & Martmez-Rogeil 
(1999)T c ff -color fitting functions, we took the four colors (B — V), 
(J — K), (V — I c ), and (V — K) as a reference for our calibra- 
tion, leading to constrain T c g for each stars in our sample within an 
estimated error better than ±100 K (see Sec. 4.1), along the whole 
spanned temperature range (3300 ^ T e g ^ 5000 K). 

The fiducial temperature allowed us to shape the unsampled 
portion of the SED at UV and far-IR wavelength by assuming a 
black-body emission independently rescaled such as to connect the 
short and long wavelength edge of the observed spectra. As shown 
in Sec. 4.2 (see also Fig.|7j, under the black-body assumption, the 
internal uncertainty in our temperature scale only impact by a few 
0.01 mag uncertainty in the inferred bolometric magnitude of our 
stars. In any case, by fully neglecting any unsampled spectral con- 
tribution, our data would be overestimating Mboi by at most 0.3 
mag. 

Making use of our new database, we have been able to draw 
a convenient set of fitting functions for the BC vs. T c ff , valid over 
the interval 3300 ^ T cff < 5000 K (see Sec. 5.1, eq.@). Similar 
relationships for BC vs. stellar colors cannot be straightforwardly 
derived (eq. O, especially for the (B — V), which shows a strong 
saturation effect for stars cooler than 3700 K, in co nsequence of the 
interv ening TiO absorption at visual wavelength ( Kucinskas et al. 
l2005h . In assessing properties of such very cool stars, however, 
one has also to consider that our sample is strongly biased against 
high-metallicity values as only the red giant branch of NGC 6791 
([Fe/H] = ±0.4) hosts stars with T off < 3700 K. 

Thanks to the wide [Fe/H] range spanned by G stars in the 
five clusters considered here, we explored the possible BC depen- 
dence on stellar metallicity. As far as atomic transitions prevail as 
the main source of metal opacity in the spectra of relatively warm 
(T e ff > 4000 K) stars, our data confirm that no evident trend of 
BC with [Fe/H] is in place (see Fig. [9}. In other words, two red gi- 
ant stars of the same effective temperature but different [Fe/H] are 
virtually indistiguishable in the values of BC to V and K bands. 
Things may be different, however, for the B (and even more for 
U) magnitudes, where the blanketing effects are more and more 
severe. In fact, Fig. QT] clearly shows that metal-poor stars display 
a "bluer" (B — V) compared to corresponding metal-rich objects 
with the same T e fj. This leads us to conclude that a drift may be 
expected for BCb such as BC'b oc —0.10 [Fe/H] among stars 
with fixed value of T c ff 

To consistently verify our calibrations, we have shown in 
Fig-EDplots of BCv and BCk vs. colors and T e g, respectively, by 
comparing with different theoretical and empirical calibrations cur- 
rently available in the literature. As far as theoretical predictions are 
concerned, it seems that the H00 models are the best ones matching 
our data in every relationship. This feature is not a surprising one, 
however, given the recognized intention of the H00 calculations to 



match M stars via "ad hoc " tuning of molecular opacity. Actually, 
this successful comparison may add a further piece of evidence, all 
the way, to the persisting limit for theory to independently assess 
the modelling of cool stars. 
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